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Abstract 

We show, using the periodic Anderson model, that the local spin self-energy approximation, as 
implemented in the extended dynamical mean field theory (EDMFT) , results in a first order phase 
transition which persists to T = 0. Around the transition, there is a finite coexistence region of 
the paramagnetic and antiferromagnetic (AFM) phases. The region is bounded by two critical 
transition lines which differ by an electron-hole bubble at the AFM ordering wave vector. 

PACS numbers: 71.27.+a, 71.10.Hf,72.15.Qm,75.20.Hr 
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I. INTRODUCTION. 



Competing Kondo and RKKY interactions in Heavy Fermion materials induce a quantum 
phase transition 111 M near which various deviations from the Landau-Fermi hquid behav- 

ri 

ior are observed experimentally j3|. Among the well-studied heavy Fermion compounds is 
CeCuQ-xAux 0| on which neutron scattering and magnetometry experiments showed 
that, in the quantum critical region, the spin susceptibilities, both the homogeneous one 
and that at the antiferromagnetic (AFM) ordering wave vector, followed, 

X-\q,T) = [T^ + e"{q)]/C (1) 

with T the temperature, 6{q) a momentum dependent function which is a measure of the 
distance from the critical wave- vector, and C the Curie constant. In the experiments it was 
found the exponent a ~ 0.75, unlike a = 1 in the standard Curie- Weiss law. This same 
behavior was found, within experimental error, to be followed by the neutron scattering data 
taken at the other wave vectors. The disentanglement of the temperature and momentum 
dependences in the inverse spin susceptibilities led to the suggestion that the self-energy 
of the spin-spin interaction be local in space and correspond to the frequency-dependent 
part of the observed The theoretical formulation of this observation turned out to be 
the extended dynamical mean field theory (EDMFT). 

The EDMFT is a method developed to study, within the local self-energy approximation, 
correlated electron systems in the existence of non-local interactions [li l^]' which, in the 
context of heavy Fermions, is the RKKY interaction. It allows the dynamical screening of 
the bare interactions. As a result, EDMFT is able to describe the competing RKKY and 
Kondo interactions in a more balanced way than the original DMFT. 



EDMFT has been app lied to study the heavy Fermions via the Kondo 



and Anderson lattice 



12 



models. Early EDMFT studies [6|, 0, [l^ approached the heavy 
Fermion quantum phase transition (QPT) by following the paramagnetic (PM) solution 
until where it ceased to exist and the spin susceptibility diverged. However, the absence 
of the AFM phase in this scenario makes it difficult to judge if the critical behavior is 
associated with a continuous transition or the spinodal point of a first order transition. 



2 



(a) 




FIG. 1: (a) Hypothetical phase diagram of the periodic Anderson model [see Eq.(j2]) below]. There 
is a continuous phase transition between the PM and AFM phases and it ends up with a QCP. 



(b) and (c 



r > 
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Two scenarios of the EDMFT phase transition, which was found to be first order at 



13l |. In the figures, the Jd line is where the AFM solution disappears. At the locus 
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121. 



there is a finite jump in the magnetization which decreases with decreasing temperature 
The Jc2 line is where the spin susceptibility at the AFM ordering wave vector diverges. In the 
region in between the two phases coexist and the first order transition is represented by a dashed 
curve. Panel (c) is a sketch of the results presented in Ref. [l^ . According to the calculation, the 
Kondo temperature Tk is near the location where the Jd and Jc2 lines become closest. The lowest 
temperature reached in Ref.jl^ is T = 0.25rR-, which is represented by a horizontal dotted line in 
(c). 



To clarify this important issue, numerical studies of the phase transition from both the 



12 1 . In the solution of the 



PM and AFM sides were carried out at finite temperatures ^ 
periodic Anderson model (PAM) 12], two different transitions were found (Jd and Jc2 
lines defined in Fig^ which bounded a region where the PM and AFM phases coexisted. 



Similar behavior was also found in the solution based on the Kondo lattice model 
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This strongly indicates a first order phase transition, at least for T > 0. 



There are important questions, though, remain unanswered. First, given the totally 
different behaviors along the mean field transition lines, it is interesting to compare 
and contrast the physical meanings of the two. Unlike at the Jc2 line where the spin 
susceptibility at the AFM ordering wave vector becomes critical, it is unclear from the 



EDMFT calculation itself 



12i | which response function is driven critical, even though 
critical slowing down was experienced. Second, there are concerns with regard to a possible 
quantum critical point (QCP) where the Jd and Jc2 lines merge ^see FigHJ^b)]. As a result, 
a novel quantum critical behavior may occur. Existing analysis [14J can not rule out such 
a possibility. Besides, although our numerical results with T > 0.25TV does not seem to 



12| may not be low 



support this scenario [see FigiU^c)], the temperatures reached in Ref. 
enough to be conclusive. The current paper is contributed to clarify these issues, which are 
all related to the local spin self-energy approximation. 

In Seciniwe introduce the EDMFT approximation on two sublattices via Baym-Kadanoff 
functional which is then used in Sec illll to formulate the instability criteria. Technical 
details of these two parts are presented in Appendices [X] and |B1 respectively. Sec JIVI 
contains conclusions and further discussions. 



II. EDMFT FORMULATION OF THE PERIODIC ANDERSON MODEL 
A. The Periodic Anderson Model 

We study the periodic Anderson model (PAM) with the local f-moments forming a hy- 
percubic lattice in dimensions: 



/ 



+ f/ E Kr - 1/2) K - V2) + E sLsi. (2) 

An RKKY interaction is introduced explicitly between the z-components of the nearest- 
neighboring f-electron spins, S^^ = nf^ — n{^^. After intergrating out the c-electrons and 



introducing a Hubbard-Stratonovic field to decouple the interactions, we obtain: 

- \ f dr f dT'Y.(^{i,T)[D,]-\cr\i'T'W,T')- f dTY,<P{i,T)Sl{i,T). (3) 
2 70 Jo ^ Jo Y 

In this action the hybridization broadened f-band is described by the free Green's function: 

Go{k, ipn) = \ipn + IJ' - Ef - ■ (4) 

V Wn + f^- eg/ 

with pn = (2n + l)7r//3. The free Boson Green's function is given by: 

J d 

Do{k, iujri) = Do{k) = -f/ H — ^ cos hi (5) 

^ i=i 

with Un = 2nn/i3. Since the bare interaction is instantaneous, the r.h.s. of Eq.(jHI) is 
frequency independent. 



B. EDMFT via Baym-Kadanoff Formulation 

We formulate the EDMFT via Baym-Kadanoff functional Qj: 

rBK[G,D,m] = Tr InG - Tr G^^G - ^Tr In D + ^Ti D^^D 

+ KnDQ^m + ^edmft[Gioc, Dioc, m]^ (6) 

G (D) is the full electron (Boson) Green's function, m = (0). The EDMFT approximated 
potential ^edmft ^ two particle irreducible (2PI) functional of the local Green's functions 
only. Since the action Q contains just a spatially local interaction vertex, ^edmft can be 
written as a summation over the local contributions. On a bipartite lattice with sublattices 
A and B, this potential is given by: 



rp 

^EDMFT[Gloc,Dioc,m] = - drY^m^J , T)aG^{jT- \jt) + Y, A[Gjj,a, Djj] + Y,'^ B[Gll,a, Du] 

° j> jeA l€B 

(7) 
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where the functionals \1/ contain second and higher order diagrams in terms of the interaction 
vertex. To solve the AFM phase with the single impurity EDMFT, we need to further assume 
Csee Ref.b 



(8) 



Here the translational invariance within each sublattice is utilized. In the PM phase the 
electron Green's functions are spin independent and the assumption is still valid. 



The Baym-Kadanoff functional gives physical solution at its stationary point. As a result, 
we have, 

/ 



AB 



ipn + fJ'- Ef 

ipn + li- Ef 



y -eg exp{-ik^) ipn + fi 



^aiiPn) 

T.^„{ip„ 



(9) 



(10) 



mn 



(11) 



Several remarks are in place. First, due to the sublattice structure, the electron Dyson 
equation Q is in a 2 x 2 matrix form. The electron self energy. 



/ . N dcf EDM FT I 
l^xAWn) = 77^ 7— ^Ijgx, 



6G,,A^Pn)"^^' ^^^^ 

with X = A,B, is local in space. Due to translational invariance, we can neglect its spatial 
coordinates. From Eq.(jHj). 



T,A,a{iPn) = T.B-a{Wn)- (13) 

As a result we are allowed to suppress the sublattice index of the self-energy in Eq.Q. 
In the PM phase, the self-energies are spin independent and the equation reduces to: 



G ^{k, ipn) = Go ^{k, ipn) - ^{ipn)- 



(14) 
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Second, in the Boson Dyson equation (fTUj) . the self-energy is defined as: 



dDjj{lUJn) 

Eq. lfTIH) carries a scalar form because the local Boson self-energy is the same on both 
sublattices due to symmetry. Finally, from Eq. ljllj) . the physical order parameter m is 
time independent and its momentum dependence, according to Eq.®, is restricted to 
Q (tt, ■ • ■ , vr) for both the AFM (m ^ 0) and the PM (m = 0) phases. In EDMFT, 
we solve the self-energies using an effective impurity model under certain self-consistent 
conditions.. [See Appendix IX] 



III. INSTABILITY CRITERIA 

A. Instability Criterion of the AFM phase (Jd hne) 

The general instability criterion against the formation or disappearance of a static spin 
density wave of wave vector Q is given by: 

-1 def d'^^BKjG, D, m] _ ^ .^g. 
dm* dm 

where m = m{Q,iO). Here the total derivatives are taken on the physical manifold of the 
Baym-Kadanoff functional defined through Eqs.(P|)-(|lip. As a result, the criterion becomes, 
(see Appendix iBl for details) 

Xt = D,\Q) -fdrj: {-ir'Uj^^^^[{r\T), (0|0)] = (17) 

•^0 I -I 

where 



a, 6=1 



'[(n|r0,(r2|r^)]= Xn^^i " XnLp + X."! + [(nkO, (r^jr^)] (18) 



In the above equation Xq q is an electron-hole bubble evaluated with the full Green's 
functions at the wave vector Q. Xo,imp is a smiliar bubble obtained via the full impurity 
Green's function. Ximp is a four point response function of the impurity model. Vo, 
which depends on T>q, is the bare interaction in the impurity model. Eq. (jl7|) gives the 
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general EDMFT instability criterion without further approximation and applies to the 
Jci line where the AFM solution at Q = (vr, ■ ■ ■ , vr) disappears. The existence of the 
electron- hole bubble at the ordering wave vector in Eq.(|18p reveals the fact that even in 
the infinite coordination limit where the mean field method becomes exact, there is still a 
non-vanishing momentum-dependent contribution in the effective spin susceptibility. The 
matrix inversions in Eq. lfTHj) involve matrices labeled by four time coordinates, two for the 
row and two for the column, respectively. As a result, this expression is generally very 
complicated and can not be further simplified. 



B. Instability Criterion of the PM phase {Jc2 line) 

In the EDMFT of the PM phase, the effective susceptibility given in Eq. fllHj) contributes 
directly to the spin self-energy and, as a result, should be local in space. This means 
we need further to restrict xo to be local. However, from the EDMFT self-consistency that 
the local Green's functions on the lattice equal to the impurity ones, Xo,ioc=Xo,imp- Hence 
the instability criterion becomes 



1 def 
A Jc2 



D,\Q)- ['dTj:i-ir''UjJir\r),m] = 
Jo 



(19) 



a,b=l 



where 



-1 



[in\r[),ir,\r^)] 



— 1 I 'T^ 

'X.imp ' 



[in\r[),iT,\r^)] 



(20) 



The special form of the matrix Vq [see Eq. fjB10|) ] allows us to carry out the matrix operations 
explicitly and obtain. 



DoiQ) = + x:'m = n-'(^o) 



(21) 



with 



def 



Xzz[r) = (5/(r)5/(0)) = [Ximp,G^G^ - Ximp^G^G^ - Ximp,GiG^ + Ximp.GxcJ [(^k), (0|0)] (22) 

The last equality in Eq.()21|) is derived by an identity 16]. Comparing with the Boson 
Dyson equation (fTU|). we see that the above instability criterion is identical to the Stoner 
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criterion in which the divergence of the magnetic susceptibihty at Q = (vr, ■ ■ ■ , vr) signals 
phase transition [l7|. 



IV. CONCLUSIONS AND DISCUSSIONS 



We have derived in this paper the phase instabihty criterion of the EDMFT solution to 
the periodic Anderson model for both the AFM and PM phases. The generic instability 
criterion (which applies to Jci line in Fig^ involves an effective spin susceptibility, Eq. (|T7j) . 
and is different from that used to determine the transition line (Jc2) bounding the PM phase, 
Eq. (fnH) . The difference is in an extra electron-hole bubble at the AFM ordering wave- vector 
in the former. This bubble is momentum dependent and survives in the infinite coordination 
limit. As a result, at the locus where one of the phases reaches the instability condition, 
the other one remains stable. This explains the phase coexistence. It persists to T = 
since the electron-hole bubble remains non-zero. This is consistent with what we obtained 
numerically in Region II of Fig^c). We should point out, though, at dimensions d > A 
and temperatures T > TKondo [Re gion I in FigHJ^c)], the difference between the Jd and 
Jc2 lines becomes negligibly small [l8|. This is due to the spatial correlation becoming 
weaker at higher dimensions and temperatures. We note in passing that no matter 
which criterion is satisfied, the divergence of the corresponding effective spin susceptibility 
at the AFM ordering wave vector naturally results in the divergence of the local spin 
susceptibility as long as the spin fluctuations are two dimensional 0, Q|. This is a re- 
sult of the dimensionality and has nothing to do with the spin self-energy being local in space. 



The true mean field transition is thus first order and lies between the Jd and Jc2 lines 
where the free energies of the two phases cross. Physically, the two sublattice EDMFT 
(as applied in the AFM phase) contains in its instability criterion an electron-hole bubble, 
which serves as a rough description of the feedback from the electron-hole excitations 
to the spin response. However, this feedback does not appear explicitly in the EDMFT 
self-consistency, which is evident from what we described in Appendix ^ As a result, 
the EDMFT spin susceptibility, which is different from the physical one in the instability 
criterion, Eq. ()17|l . does not experience any singularity as the Jd line is crossed. On the 



9 



other hand, the homogeneous EDMFT (as apphed in the PM phase), contains the same 
singular behavior in the spin response as that in the instabihty criterion, Eq. (pi?|) . As a 
result, when the phase boundary is approached, EDMFT is able to adjust self-consistently 
to reflect the singular behavior in the spin channel. However, as we have already noted, 
the problem on this side is that the feedback from the non-local electron-hole excitations 
is totally missing. So both the transition lines contain unphysical features, and neither 
of them, as far as the critical properties are concerned, is close to the true transition. A 
related issue, which concerns the critical exponent a in Eq.(^ along the Jc2 line, further 
supports our conclusion. It was shown that at T = on the Jc2 line, the critical frequency 
dependence could not develop a sublinear form [iQ^ . 



After all, it is not a surprise that, although it works well qualitatively in describing many 
other physical properties 3|; the EDMFT fails to capture the right phase transition. This 
is certainly one of the issues one needs to improve over the mean field approach. Given 
what we have concluded in this paper, it seems important that one needs to find a way 
allowing proper feedback from the electron-hole excitations, which is spatially non-local, to 
the f-electron spin response. A natural way to proceed is to combine the EDMFT scheme 
with the random phase approximation (RPA) |2Q|]. In this combination, the spin self-energy 
contains the local EDMFT part together with the non-local RPA part. This is a desirable 
feature as one can see from the EDMFT instabilit y c riterion Eq. ()18|l . Besides, the scheme 
is derivable from the Baym-Kadanoff functional [2i|. Of course, with the new scheme, 
the instability criterion itself is modified and its implication to the heavy Fermion phase 



transition has not yet been explored. A different route is to utilize the cellular DMFT 



2l|. 



To this end, a two impurity Anderson model subject to the DMFT self-consistent electron 



221. In this formalism, the RKKY interaction 



bath results in a qualitative improvement 
is generated dynamically, instead of being added in by hand as in Eq.(l2)). The spin 
susceptibility across the two impurity sites, which contains the corresponding electron-hole 
bubble as the lead order contribution, renders a limited momentum dependence and turns 
out to be essential to the improvement. 
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APPENDIX A: EFFECTIVE IMPURITY MODEL AND EDMFT SELF- 
CONSISTENCY 

To obtain the local self-energies, we need to solve an effective impurity model: 

Jo Jo „ 

-\f dT rrfr'0(O,r)Po-i(r-r')0(O,r')- dT<j>{0,T)S[{0,T). (Al) 
/ Jo Jo Jo 

The mean field Weiss functions Qoa and Vq are decided by the following self-consistent 
conditions: 

So.\Wn) = E GAk, ^Pn)]-' + (A2) 

it 

T^o\^^n) = E i^n)V + (A3) 

k 

with Go-{k,ipn) and D„{k,iUn) given by Eqs.© and (fTUj). respectively. The self-energies 
are, 

Sr(^Pn) = [<3o.r\Wn) - [GTT\Wn) (A4) 
Yi'^^ilUOn) = m-\lUJn) - [D''^T\^^n) (A5) 

where the impurity Green's functions G^^^ and D*"^? are obtained by solving the effective 
action ()A1|) . In Eqs.Q and (fTUj) . we need to use the lattice self-energies which are usually 
assumed to be the same as the impurity ones in the disordered phase. In the ordered phase, 
the electron self-energy on the lattice is different from that of the impurity model by a 
Hartree term, while the Boson self-energy is still the same. 
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j:^itPn) = T^r^Wn) - a[Vom - (f/ - JRKKY)]{Si) (A6) 

Yi{iun) = (A7) 

The meaning of Eq. ()A6|l is that we need to replace the Hartree self-energy of the impurity 
model by that on the lattice, using the electron magnetization. This procedure is related 
to Eq. ()ll|) . which would otherwise introduce a third self-consistent equation. With this, we 
have presented a complete self-consistent loop. 



APPENDIX B: DERIVATION OF THE INSTABILITY CRITERION FOR THE 
AFM PHASE 

We derive here the instability criterion specific to the periodic Anderson model. From 
the general condition, Eq. lfTHIl . together with Eqs.(P)- (fTT|l . we obtain. 



dx 



dm* dm 



BK 



BK 



dm* dGa^{x)dG^y{y) 



dm* dGa^{x)dm 

f , dD(x) d'^T 
+ / dx- 



+ dx 



dD{x) d^TsK 



dm* dD{x)dm 

BK 



+ / dx 



dm* dD{x)dGay{y) dm*dGcryiy) 

BK 



dD{x) d'^r 

BK , (J i- } 



(Bl) 
(B2) 
0. (B3) 



dm* dGfj^{x)dD{y) ' J dm* dD{x)dD{y) ' dm*dD{y) 

We used x = {Rj,T\Rji,T'), (similar for y) and J dx = lo dr Jq dr'. Summation over 
the repeated spin indices is implied. Solving dG/dm* and dD/dm* from Eqs. ()B2|) and ()B3|1 . 
and substituting them in Eq. ()Bl|) . we obtain: 

d^TBK[G,D,m] 



dm* dm 



dx / dy d'TBK/dm*dG^^{x),d'TBK/dm*dD{x) 



d^TsK/dG^^dG^^ d^TsK/dG^^dD 
d^TBK/dDdG^y d^TBK/dDdD 



-1 



{x,y) 



d^TBK/dG^yiy)dm 
d'^TBK/dD{y)dm 



0. 



(B4) 



Using Eqs.® and ((Tj), we find d'^T bk / dG^{jTi\j'T[)dm{i, t) = —a5ij6iji6{T — ti)5{t — r{) 
and d'^TBK/dD{jTi\j'T[)dm{t,T) = 0. Besides, d^TBK[G, D,m]/dm*dm = Dq\Q). So we 
have: 



12 



(3D^\Q) = ^jY.fdnY. f dT2 exp{-tQ ■ i?, J exp(iQ ■ 4) 
N Jo Jo 



x(l,-l,0) 



p(2) pv^; p^ 

G-^D 



.(2) 



p(2) p{2) -p^ 

CjGi GiGi G|D 



J2 

.(2) 
^(2) 



.(2) 



.(2) 
DG| 



.(2) 



-1 



[(jinljin), (j2T-2b2T-2)] 



-1 

V 0/ 



(B5) 



with r5,'l.[(jiri|jlrO, (J2r2|j^r^)] = 52rBi^/9X(jiri|j(r09r (j2r2|j>^) for X, F = To 
solve the matrix T^'^\ we use again the Baym-Kadanoff functional ® and obtain: 



(2) 



(B6) 



where 



XO,XY [{jlTl \j[r[), {]2T2 I JsTa)] 



dof 



/^(jin|j2r2)/^(jir{|j^r^)+Z}(jiri|j^r^)D(jlr{|j2r2), X, F = 

0, else 

(B7) 



<^>?l-[(jin|jirO,(j2r2|j>D] 



def 



52$ 



EDMFT 



dX{j,T,\j[T'^dY{j2T2\3'2r^) 



(B8) 



(2) 

$)(^y, same as ^edmft, contains only propagators local in space and is 2PI in separating 



the external legs labeled by 1 and 1' from those 2 and 2'. It follows then 



ft'^l l{hn\ATi},U2r2\A<)] 

L IXY 



-xoLp + XiL + ^0 [(nkO, (^-2^2)] 



XY 



(B9) 

Here Xo,imp is similar to that defined in Eq. ()B7|) except being local in space. We also defined: 



^o[(n|rO, {r2\T2)] = 6{ti - t[)5{t2 - Tg) 



Voin~T2) ~Voin-T2) 
-Voin-T2) Vq{ti-T2) 





(BIO) 
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Xirnp,xy[(0ri|0r0,(0r2|0r^)] (T. : 0^(0, n) 0^(0, r^) :: 61(0, 61(0, r^) ■) (Bll) 
where Ox = Cfj (0) if X = Go- The instabihty criterion becomes: 



Dq\Q) = J dTexp{-iQ ■ Rj^)exp{iQ ■ Rj^) 



a,b=l 



(B12) 



where all the four terms in the square parenthesis are 3x3 matrices and after matrix 
inversion, only the first 2x2 block contributes. It should be noted that the matrices 
are also labeled by the two pairs of the space-time coordinates and any matrix operation 
should take these into account. As a result, e.g., in Eq. ()B12|) the full matrix, labeled by 
lUi'^ilji'^i)^ U2T2\j2T2)]^ should bc inverted first and only after that, we set the labels to be 
[(Jir|jir),(j20|j20)]. 

Finally, since in Eq. ()B12|) . Xo is the only term contains spatially non-local contributions, 
the Fourier transform over the lattice coordinate can be taken into the matrix inversion, 
which gives: 



(B13) 



where 



-1 



imp 



[in\T[),{T2\T^)] 



(B14) 



This gives an instability criterion consistent with the EDMFT Baym-Kadanoff functional, 
Eqs.dni) and ((Tj), without any further approximation. 



As it turns out, we need further to assume xo be spatially local in Eq. ()B12|) . in order to 
describe the PM phase. In such a case, Xo Xo./oc, we have (1) the momentum dependent 
phase factors in Eq. ()B12j) cancel out and (2) Xo,imp cancels Xo,ioc due to the EDMFT 
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self-consistency that the local lattice Green's functions equal to the impurity ones. 
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